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Abstract. Recent data from high statistics experiments that have measured the modulus of the pion 
electromagnetic form factor from threshold to relatively high energies are used as input in a suitable 
mathematical framework of analytic continuation to find stringent constraints on the shape parameters of 
the form factor at t = 0. The method uses also as input a precise description of the phase of the form factor 
in the elastic region based on Fermi- Watson theorem and the analysis of the nir scattering amplitude with 
dispersive Roy equations, and some information on the spacelike region coming from recent high precision 
experiments. Our analysis confirms the inconsistencies of several data on the modulus, especially from low 
energies, with analyticity and the input phase, noted in our earlier work. Using the data on the modulus in 
a region of stability, taken as (0.65, 0.70) GeV, we obtain, with no specific parametrization, the prediction 
(rj) G (0.42, 0.44) fm 2 for the charge radius. The same formalism leads also to very narrow allowed ranges 
for the next higher-order shape parameters at t = 0, with a strong correlation among them. 
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1 Introduction 

The pion electromagnetic form factor is an important quan- 
tity that encodes information on the structure of the strong 
interactions. At high spacelike momenta, it acts as an ex- 
cellent testing ground for studying the onset of pertur- 
bative QCD for exclusive quantities, while at low ener- 
gies it is an important laboratory for the study of chiral 
symmetry breaking. It plays a major role in the precision 
tests of the Standard Model, through its contribution to 
the anomalous magnetic moment of the muon, which is 
one of the most precisely measured observables in particle 
physics. 

The derivatives of the form factor at t = present in 
the Taylor series expansion 



F(t) = 1 + -(r*)t + ct 2 + dt 3 + 



(1) 



are quantities of interest for testing the expansions of Chi- 
ral Perturbation Theory (ChPT) and the lattice calcula- 
tions of the pion form factor. These shape parameters, 
especially the charge radius squared (r 2 ), are the object 
of many investigations. However, the origin is not directly 
accessible to experiment. Measurements near the origin 
EL Hi & EL Hi m the spacelike region t < have been 
exploited in the past to yield the charge radius and the 



higher shape parameters. Recent high statistics experi- 
ments @, 0, H, d, [H| have measured the modulus in 
the timclike region on the unitarity cut. Since this exper- 
imental information is available only quite far away from 
the origin, new and suitable techniques of extrapolation 
have to be invented in order to obtain a reliable extrac- 
tion of the shape parameters from these data. 

The analyticity properties of the form factor are very 
useful in performing this extrapolation. From causality 
and unitarity it is known that F(t) is an analytic func- 
tion in the complex i-plane, with a cut starting at the 
unitarity threshold t + = AM 2 and running to infinity. 
The Fermi- Watson theorem relates the phase of the form 
factor in the elastic region, t+ < t < t[ n , where t- m is the 
first inelastic threshold, to the P-wave phase shift of 7T7T 
scattering, which is now known to high precision due to 
studies based on Roy equations for the 7T7t clastic ampli- 
tude BGl El Ellil. 

Analyticity was exploited in the past by means of var- 
ious types of dispersive representations, which express the 
function in terms of either the imaginary part, the phase 
or the modulus on the cut. However, none of the standard 
dispersion relations has complete input, and their appli- 
cation requires ad-hoc assumptions. Specific parametriza- 
tions of the data were also frequently used, but they are 
affected by the known "instability" of analytic continu- 
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ation [7j} if they are extrapolated outside their original 
range of validity. 

Recently, more sophisticated methods of complex anal- 
ysis were applied in a formalism known as "theory of uni- 
tarity bounds" (HQS. The 

main feature of this approach 
is that it allows the optimal implementation of the phase 
and modulus information available on the unitarity cut. 
Using mathematical techniques belonging to the analytic 
interpolation theory on Hardy spaces, one can include as 
input also values of the function or its derivatives at dis- 
crete points inside the analyticity domain. This input is 
shown to lead to strong constrains on other values of the 
form factor, like the shape parameters at the origin or the 
values on the spacelike axis. 

The inherent limitation of this approach is that it pro- 
duces only bounds on the quantities of interest, rather 
than definite determinations. Moreover, when the experi- 
mental errors are large, there is no significant advantage in 
adding more and more input information, as the bounds 
do not improve. On the other hand, due to the increased 
accuracy of the input observed recently, the bounds can be 
very stringent, competing in precision with experimental 
data or theoretical predictions. 

Applications of this formalism to the pion form factor 
were reported in [l8j] and more recently in 0, Hi], [HJ ■ In 
[20( | . the phase in the elastic region t < ti n , where t- m is 
the first important inelastic threshold set by the produc- 
tion of an ujTT pair, together with the modulus above t- m 
used in an averaged way and some spacelike information, 
were exploited for deriving constraints on the higher shape 
parameters c and d appearing in the expansion (TT|). 

In a more recent work [22j , the same input was used for 
deriving bounds on the modulus in the elastic region, be- 
low the inelastic threshold. The results provided nontrivial 
consistency checks on the recent experimental data @-[Hl 
on the modulus from e + e - annihilation and r-decay ex- 
periments. In particular, at low energies the calculated 
bounds offered a more precise description of the modulus 
than the experimental data. 

In [2(| HH, HH j the experimental data on the modulus 
below the inelastic threshold were not included as input. 
The results showed that the knowledge of the phase in 
the elastic region has strong implications for constraining 
the shape of the form factor below and above the uni- 
tarity threshold. The inclusion of the experimental mea- 
surements on the modulus of the form factor below the 
inelastic threshold is the object of the present work. More 
precisely, we ask what is the import of these measure- 
ments on the radius (r 2 ) and the higher derivatives c and 
d at t = 0. Thus, the present work extends our previous 
analysis of the shape parameters in (20l |. where the infor- 
mation on the modulus below tj n was not used. It offers a 
novel way in which the analyticity methods can be used 
to complement the work in Refs. [1(1 HJL H3 and provides 
a further consistency check on the experimental data. 

There is a rich discussion in the literature on the de- 
termination of the shape parameters. On the theoreti- 
cal side, a fit based on ChPT to two-loop accuracy for 
r decays gives (r£) = (0.431 ± 0.020 ± 0.016) fm 2 and 



c = (3.2 ± 0.5 CX p ± 0.9 th cor) GeV -4 [13 . The pion form 
factor has been calculated in two-loop ChPT which gives 

(li) = (P4 37 ± °- 016 ) f m 2 and c = ( 3 - 85 ± °- 60 ) GeV~ 4 
[24|. In [25|. the values (r 2 ) = (0.452 ± 0.013)fm 2 and c = 

(4.49±0.28) GeV -4 were obtained within the three flavour 
framework and at next-to-ncxt-to-leading order in ChPT. 
Recent lattice calculations with chiral extrapolation based 
on two-loop ChPT give (r 2 ) = 0.409(23)(37) fm 2 and c = 
3.22(17)(36) GeV" 4 [1|. Phase (Omnes-type) representa- 
tions with various parametrizations of the phase along 
the whole unitarity give (r 2 ) = (0.432 ± 0.001) fm 2 and 
c = (3.84 ± 0.02) GeV -4 [27|. Recently, a fit of space- 
like data with Pade approximants [28l ] gave the values 
(r 2 ) = (0.445±0.002±0.007)fm 2 and c = (3.30±0.03 sta t± 
0.33 syst ) GeV- 4 . 

Other values are (r 2 ) = (0.427 ± 0.010) fm 2 , quoted 
in Q, and (r 2 )^ = (0.711 ± 0.009) fm given in Q, both 
being obtained by fits of the data with a simple pole. The 
curvature c has also been determined from fits of experi- 
mental data with specific analytic parametrizations of the 
form factor. The value c = (3.90 ± 0.10) GeV -4 has been 
obtained in [1^] by a standard dispersion relation. A fit 
of the ALEPH data |3U on the hadronic r decay rate 
with a Gounaris-Sakurai formula for the form factor [3l| 
gives c = (3.2 ± 1.0) GeV -4 . The analysis based on the 
technique of unitarity bounds [20l ] mentioned above led to 
3.75 GeV -4 < c < 3.98 GeV" 4 . 

The next shape parameter d is much less known. The- 
oretical results from ChPT and lattice calculations are not 
yet available. The value d = (9.70 ±0.40) GeV -6 has been 
obtained from fits of data by means of usual dispersion 
relations [1^] , while the Taylor expansion of the Gounaris- 
Sakurai parametrization [3(|, leads to d — 9.80 GeV -6 . In 
we derived the range 9.91 GeV -4 < d < 10.46 GeV -4 . 

Except for the ranges predicted in [2QJ, the previous 
determinations of the shape parameters came from spe- 
cific parametrizations, like polynomial expansions, Pade 
approximants of the Gounaris-Sakurai model for the dom- 
inant p-pole. In contrast, the predictions of the present 
work are parametrization-free in the sense that we do not 
rely on specific analytic expressions of the form factor or 
its modulus. The determination follows from general prin- 
ciples of analiticity with the input being the phase in the 
elastic region and conservative modulus information on 
the unitarity cut. The method also avoids the instability 
problems inherent in the analytic continuation of specific 
parametrizations . 

The scheme of the paper is as follows: in Sec. [3J we 
present very briefly the mathematical formalism whereas 
in Sec.[3]we describe the information used as input to our 
work. In Sec. [H we present the results of our investigations. 
In subsection 14.11 we present the analysis of the charge 
radius (r 2 ) and in subsection [472] the analysis of the higher 
shape parameters c and d. In Sec. [5] we summarize our 
results and present some conclusions. 
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2 Method 

We use the Fermi- Watson theorem, which states that 

Arg[F(f + ie)] = Sj(t), t+ < t < t in , (2) 

where 5{ (t) is the phase-shift of the P-wavc of 7T7t clastic 
scattering and f; n the first inelastic threshold. In addition, 
we exploit the recent experimental data on the modulus 
above f; n and the 1/f asymptotic decrease of the form 
factor predicted by perturbative QCD, by adopting the 
relation 

- / dtp(t)\F(t)\ 2 = I, (3) 

n J tin 

p(t) is a suitable positive-definite weight allowing a reliable 
calculation of the value of /. As in [20, H3, [22j we consider 
weights of the form 



Pit) 



t" 



(4) 



where Q 2 > and 6, c are taken in the range b < c < b + 2. 

The conditions ^) and ([3]) can be written in a form 
that allows the application of the mathematical interpola- 
tion theory for analytic functions [H, [33[ . We first exploit 
([2| by introducing the Omncs function 



0(f) = cxp 




dt'- , 



(5) 



where 5(t) = 6{(t) for t < t- m , and is an arbitrary function, 
sufficiently smooth (i.e., Lipschitz continuous) for t > U n . 
The crucial observation is that the function h(t), defined 
by 

F(t) = 0{t)h{t), (6) 

is real for t < t; n , i.e. it is analytic in the t-planc cut only 
along t > t m . Furthermore, from it follows that h(t) 
satisfies the condition 



dtp(t)\<D(t)\ 2 \h(t)f 



(7) 



This relation can be written in a canonical form by per- 
forming the conformal transformation 



2(f) = 



ti n H~ yt>in t 



(8) 



which maps the upper (lower) lip of the branch-cut [fj n , oo] 
to the upper (lower) half of the unit circle in the complex 
z-plane, and the cut f-plane onto the interior of the unit 
circle \z\ < 1, the real line [— oo, 0] to [—1,0] and [0,fi n ] 
to [0, 1]. We introduce then two outer functions, i.e. func- 
tions analytic and without zeros in the unit disk \z\ < 1, 
defined in terms of their modulus on the boundary, re- 
lated to y/p(t) |df/d2(f)| and \0(t)\, respectively 0Q1]. 
In particular, for weight functions of the form the first 



outer function w(z) can be written in an analytic closed 
form in the z-variable as UM 



W(z) = (2y/t~) 



1+b-c 



(1-z) 1 / 2 (l + z(-Q 2 )) c 
(1 + z) 3 / 2 - c + b (1 - zz(-Q 2 )) c ' 



(9) 



For the second outer function, denoted as uj(z), we use an 
integral representation in terms of its modulus on the cut 
t > ti n , which can be written as [TH, [l^| 



uj{z) = exp 



\Ain " i{ Z ) 



ln|C(f')|df' 



h in \/F^f~(f - t{z)) ) ' 
( 10 ) 

where t{z) is the inverse of z = z(t), for z(t) defined in 



Further, we define a function g(z) by 

g(z) = w(z)u;(z)F(t(z))[0(t(z))}- 1 . 



(11) 



Then ([7]) can be written in terms of the integral of \g(z) 
on the boundary (z = e 16 ) as 

r-2n 



1 



d0\g(e«)\' 



I. 



(12) 



The L 2 -norm condition ([T2l leads to what is known as the 
Meiman interpolation problem 32] . It consists of finding 
the most general rigorous correlations between the values 
of the function and its derivatives inside the unit disk \z\ < 
1 consistent with the relation (|12p (for a proof and older 
references see [Hj]). 

For instance, with techniques of complex analysis, one 
can show that (fT2"|) implies the following determinantal 
inequality : 



I 

6 



6 

z l 


6 • 

[ziz 2 ) K 


■ £n 

(ziz n ) k 


l-Z 2 

(ziz 2 ) k 


1 - z x z 2 
(z 2 ) 2K 


1 - ZiZ N 

(z 2 z N ) K 


1 - ZxZ 2 


l-Z 2 


1 - z 2 z N 



N 



(z 1 z N ) K (z 2 z N ) 



K 



■2K 
N 



1 - ZlZ N 1 - Z 2 Z N 

where the auxiliary quantities 

A'-l 



> o, 



(13) 



-N 



K-l 



I = I ~ X] 9k, £n = g{z n ) - 9 k2 



(14) 



fc=0 



k=0 



are defined in terms of the values 
1 d k g(z) 



k\ dz k 



= 9k, 0<k<K-l, 



g{z n ) = &i, 1 < n < N. 



(15) 



For simplicity we considered N real points z n and K — 1 
derivatives at z = 0. The details of the derivation are 
reviewed in Ref. 1191 . 
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By using (fTTjl one can express the inequality (fT3"]l as 
a quadratic constraint on the values of the form factor 
F(t) and its derivatives at specific points. From this, by 
solving simple quadratic equations, one can derive upper 
and lower bounds on the first K — 1 derivatives at t = 0, in 
terms of the other values included in the set. In particular, 
if the point t n is situated on the elastic part of the cut, 
i.e. t+ < t n < t- m , the relation (|16p writes as 

g(z n ) = w{z n )uj(z n ) \F{t n )\/\0{t n )\, (16) 

where z n = z(t n ) and the modulus |0(f)| of the Omncs 
function is obtained from ([5]) via the Principal Value (PV) 
Cauchy integral 

\0{t)\ = exp ( Ap V ^ dt'-^P-r) . (17) 
\ir J 4M 2 t'{t'-t)J 



3 Input 

We recall that the first inelastic threshold t- ln for the pion 
form factor is due to the opening of the wk channel which 
corresponds to y/tfa = + M„ = 0.917 GeV. Below t[ n 
Fermi- Watson theorem ^ relates the phase of the form 
factor to the P-wave phase shift S{ (t) determined recently 
with high precision from Roy equations satisfied by the 
7T7T clastic amplitude [l2|, EH, Q H, El] ■ We calculate the 
Onmes function ([5]) using as in put for t < t- m the phase 
shift 5{(t) from Refs. |15[ and [lg|, which we denote as 
Madrid and Bern phase, respectively. Above t- m we use 
a continuous function S(t), which approaches asymptoti- 
cally 7r. As shown in [l9j, if this function is Lipschitz con- 
tinuous, the dependence on S(t) of the functions 0(t) and 
ui(z), defined in ([S]) and (|TU|) . respectively, exactly com- 
pensate each other, leading to results fully independent 
of the unknown phase in the inelastic region. This is one 
of the important strengths of the method applied in this 
work. 

In our analysis we use the condition ([3]) with weights 
of the type ([!]). We calculate the integral defined in ((3J) 
using the BaBar data @ from t- m up to \fi = 3 GeV, 
continued with a constant value for the modulus in the 
range 3 GeV < \ft < 20 GeV, smoothly connected with a 
1/t decrease above 20 GeV. As discussed in [22j , the choice 
of the weight should lead to an accurate value for the 
integral I, providing at the same time a strong constraint 
on the high energy behaviour of the form factor. We have 
tested a large class of expressions, choosing finally in our 
analysis the weight p(t) = 1/t, for which the value of I is 

7 = 0.578 ±0.022. (18) 

It turns out that the results for the shape parameters ob- 
tained with various weights are very similar. We illustrate 
this by reporting also several results obtained with the 
weight p(t) = \j\ft. for which [22j : 

I = 0.687 ±0.028. (19) 



The uncertainties quoted in (|18|) and (fTTJl) are due to the 
BaBar experimental errors. As shown in jl8l[l9| . the bounds 
depend in a monotonous way on the value of the quan- 
tity I, becoming stronger /weaker when this value is de- 
creased/increased. In the applications, we have used as 
input the central value of / increased by the error, which 
leads to the most conservative bounds. 

We use as input also the values at several real points 
inside the analyticity domain \z\ < 1, which translate as 
t < tj n in the t-plane. Specifically, we implement the nor- 
malization condition 

F(0) = 1, (20) 

which follows from ([I} , and one of the spacelike data taken 
from [1(1]: 

F(-1.60 GeV 2 ) = 0.243 ± 0.012±g;gJ|, 
F(-2.45 GeV 2 ) = 0.167 ± 0.010+^Jf. (21) 

In addition, we have used the value of the modulus F(t n ) 
at an energy below the W7r inelastic threshold, 

\F(t n )\ =F n ± e„, t+ < t n < ti„. (22) 

with the central value F n and the error e„ taken from 
one of the recent experiments @j-[ll|. We have taken into 
account of the isospin violating effect due to the p — u> 
mixing, by applying to the data from e + e~ annihilation a 
suitable correction factor [34|, HH that we have also used in 
our previous work [22| . Therefore our results on the shape 
parameters will correspond to the exact isospin limit. 

An important remark is that, except for the normal- 
ization condition (|20[) which is exact, all the inputs that 
we use are known with some uncertainty. The treatment 
of the error channels in interpolation theory is a nontriv- 
ial problem, which requires some caution. Following the 
discussion given in [221 ] . in our analysis we have varied 
all the inputs simultaneously within their error intervals, 
taking the most conservative bounds on the quantity of 
interest consistent with the input, i.e. the largest upper 
bound and the smallest lower bound from the set of values 
obtained with specific inputs. In other words, we take the 
union of the allowed domains, obtained with specific input 
values inside the error intervals, for the quantity of inter- 
est. The procedure has been carried out efficiently with a 
combination of Mathematica and Fortran programs. 



4 Extraction of results 

4.1 (r'l) analysis 

We study first the impact of the conditions ©, ([3]), ([211)) . 
(f2"Tj) and (|2"2"j) used as input in our formalism on the charge 
radius (r 2 .). To this end we calculate upper and lower 
bounds on the first derivative appearing in the expan- 
sion ([T]), by applying the general inequality (|13[) in the 
particular case K = 2 and N = 2. More precisely, we 
use as input, besides the normalization (|20[) . one of the 
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spacelikc values (|21|) and a single experimental modulus 
|F(i ra )| at an energy squared t„. below the ujtt threshold, 
from the sets reported in [f|-[llj]. The errors on the phase, 
the spacelike value and the modulus were taken into ac- 
count as explained at the end of Sec. [3] The results are 
obtained using in © the weight p{t) = 1/t, with the value 
of / taken at the upper limit of (|18l) . 



In Fig. [TJ we present upper and lower bounds on (r^) 
as functions of the energy \ft corresponding to the input 
modulus. We have used input data from BaBar experi- 
ment @ along the whole elastic region from the threshold 
4M^ to the first inelastic threshold ij n . For convenience 
we connected the individual points by continuous lines Q 

^From Fig. [Tj one can see that the curves shown are 
not smooth: the bounds obtained for (rj) vary quite dras- 
tically from point to point, which means that they de- 
pend strongly on the experimental modulus used as in- 
put. We obtain similar fluctuations also with the modulus 
data measured by KLOE, CMD-2 and Belle experiments. 
In addition, there are several experimental data for which 
we do not even obtain real bounds for (rj): the quadratic 
equations derived from the inequality (|13[) have complex 
solutions for all the input values inside the error bars. In 
these cases (not shown in Fig. [Tj the input modulus does 
not satisfy the analyticity constraints with the phase and 
the spacelikc values used as input. Such inconsistencies in 
the data on modulus were noted already in our previous 
analysis [22j ]. 



0.46 




0.38 z_ 
0.3 



0.4 



0.5 



Bern phase [16] 
BaBar modulus [6] 
F(-1 .60 GeV 2 ) [4,5] 
J , I i L_ 



0.6 

t 1/2 [GeV] 



0.7 



Oi 



0.9 



Fig. 1. Upper and lower bounds on (r^) using as input one 
modulus value below the ujtt inelastic threshold from BaBar 
experiment [f|, as functions of the energy \fi where the mod- 
ulus was implemented. We used the Bern phase shift from [l(J 
and the spacelike value F{— 1.60 GeV 2 ) given in <|2ip . 



The upper and lower bounds shown in Fig. [T] define, 
for each input modulus at a fixed energy, an allowed in- 
terval for the charge radius (r^), which turns out to be 



1 The number of experimental points below the ujtt threshold 
reported by the BaBar experiment 0] is 221, while KLOE 7, 8] 
measured the modulus at 121 energies below ti n , CMD-2 [9l.ll0| 
at 34 and Belle [ll|] at 15 energies. 



quite narrow for some particular energies. In principle, if 
the data were consistent among each other, the final al- 
lowed domain for (rj) would be given by the intersection 
of the allowed ranges obtained with particular inputs at 
fixed energies. The range defined by this intersection will 
be limited by the smallest upper bound and the largest 
lower bound shown in the figure. Of course, if the smallest 
upper bound turns out to be smaller than the largest lower 
bound the intersection is empty. In practice this turns out 
to be the case, if we consider all the points shown in Fig. 
[TJ as expressed in Table [TJ This is seen in the left part 
of Table [2 where (r^) max is the lowest upper bound and 
(^)min is the largest lower bound upon all the experi- 
mental points below the ujtt threshold. The correspond- 
ing values for the four experiments BaBar, KLOE, CMD- 
2 and Belle are indicated separately (the data labelled 
as CLEO will be explained below). For all experiments 
( r J)max < (^)min, which clearly indicates that there are 
inconsistencies between the measurements of the modulus 
at different energies, if we consider all the data between 
the elastic threshold t + and the inelastic threshold i; n . 
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Fig. 2. Upper and lower bounds on (r 2 .) with various sets 
of input modulus data in the region of stability (0.65 - 0.70) 
GeV. We have used the Bern phase [l|J and the spacelike value 
F(-1.60GeV 2 ) from J2TJ. 



Taking into account this picture, in order to get a rea- 
sonable bound on (r^) we look for an energy region where 
the bounds are stable with respect to the variation of the 
input energy and do not exhibit big fluctuations. As can 
be inferred from Fig. [Tj and Table [TJ for BaBar data the 
region of stability can be taken to be (0.65 - 0.70) GeV. 
The same region can be taken also for the other data sets, 
KLOE, CMD-2 and Belle. The upper and lower bounds 
on (r^) obtained with various data sets in this region of 
stability as functions of the input energy are compared in 
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Table 1. Intersections of the ranges of (r 2 ) obtained with various inputs. 



Spacelike data 


\F(t)\ 


Bern phase [16J 


Madrid phase [l^j 


All points included 


(0.65 - 0.70) GeV 


All points included 


(0.65 - 0.70) GeV 


( r 7r } min 


(^7r)max 


} mill 


(^Tr)max 


(?V)miii 


(^7r)max 


( r 7r)min 




F( 1.60GeV 2 ) 


Belle 


0.4229 


0.4028 


0.4200 


0.4428 


0.4362 


0.4254 


0.4294 


0.4463 


BaBar 


0.4562 


0.4299 


0.4210 


0.4404 


0.4455 


0.4343 


0.4302 


0.4435 


CMD2 


0.4302 


0.4278 


0.4125 


0.4373 


0.4357 


0.4288 


0.4190 


0.4406 


KLOE 


0.4264 


0.4255 


0.4158 


0.4362 


0.4430 


0.4302 


0.4248 


0.4385 


CLEO 


0.4139 


0.4107 


0.4137 


0.4378 


0.4205 


0.4256 


0.4196 


0.4390 


F(-2.45GeV 2 ) 


Belle 


0.4227 


0.4027 


0.4194 


0.4418 


0.4368 


0.4250 


0.4301 


0.4455 


BaBar 


0.4562 


0.4279 


0.4206 


0.4391 


0.4455 


0.4325 


0.4309 


0.4425 


CMD2 


0.4310 


0.4257 


0.4109 


0.4360 


0.4357 


0.4285 


0.4188 


0.4395 


KLOE 


0.4267 


0.4234 


0.4148 


0.4348 


0.4430 


0.4282 


0.4250 


0.4373 


CLEO 


0.4129 


0.4104 


0.4125 


0.4366 


0.4203 


0.4251 


0.4193 


0.4379 


Table 2. Weighted average of upper and lower bounds for (r 2 ) obtained with various inputs. 


Spacelike data 


\F(t)\ 


Bern phase \16j 


Madrid phase \15j 


All points included 


(0.65 - 0.70) GeV 


All points included 


(0.65 - 0.70) GeV 






{ t tt) mill 


(^Tr)max 




(^Tr)max 


( ?I 7r)min 




F( 1.60GeV 2 ) 


Belle 


0.4152 


0.4443 


0.4209 


0.4456 


0.4187 


0.4455 


0.4269 


0.4475 


BaBar 


0.4168 


0.4467 


0.4194 


0.4439 


0.4209 


0.4474 


0.4255 


0.4458 


CMD-2 


0.4127 


0.4443 


0.4133 


0.4408 


0.4122 


0.4439 


0.4179 


0.4421 


KLOE 


0.4051 


0.4470 


0.4151 


0.4408 


0.4055 


0.4453 


0.4205 


0.4421 


CLEO 


0.4181 


0.4260 


0.4161 


0.4408 


0.4183 


0.4267 


0.4190 


0.4399 


F(-2.45GeV 2 ) 


Belle 


0.4137 


0.4433 


0.4204 


0.4448 


0.4180 


0.4446 


0.4275 


0.4468 


BaBar 


0.4155 


0.4460 


0.4187 


0.4430 


0.4203 


0.4466 


0.4259 


0.4450 


CMD2 


0.4107 


0.4432 


0.4118 


0.4396 


0.4102 


0.4428 


0.4177 


0.4411 


KLOE 


0.4013 


0.4461 


0.4139 


0.4397 


0.4021 


0.4442 


0.4204 


0.4411 


CLEO 


0.4179 


0.4254 


0.4150 


0.4396 


0.4181 


0.4262 


0.4187 


0.4387 
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Fig. 3. As in Fig. [2] but with the input value F(-2.45 GeV 2 ) 
from (pTjl . 
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Fig. 4. As in Fig. [2] using the input phase from [IB] ]. 



Figs. [5] -[5] In these figures, we show the behaviour of the 
bounds for different choices of the phase and the space- 
like datum. We observe that the bounds obtained with 



the Madrid phase [15[ are slightly shifted upwards com- 
pared to those obtained with the Bern phase [l|| . On the 
other hand, the choice of the spacelike data does not have 
a drastic influence on the bounds. We also see that the 
bounds from BaBar [|| and Belle [TTJ] are shifted upwards 



compared to the bounds obtained with the data of CMD-2 
and KLOE. 

If we restrict to the region of stability, we obtain a 
nonzero intersection of the allowed ranges for the charge 
radius obtained with individual data on the modulus. This 
is illustrated in the right part of Table [TJ where, for each 
data set, (rj) max is the lowest upper bound and (rj) m i n 
is the largest lower bound when the energy of the input 
modulus is varied along the stable region. From the values 
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listed in the Table for all data sets we infer that the final 
range obtained by the intersection of the particular ranges 
in the stable region is given by 



' 7r/min 



:0.42fm 2 , (r 2 ) max «0.44fm 2 . (23) 



<rj>« 



y ■ 



1/4 (24) 



where the sum runs over the individual energy points t n . 
The results are presented for all data sets in Tabled where 
we give separately the averages for all the experimental 
points and the stable energy region defined above. We find 
from our analysis that the effective bounds obtained from 
the weighted average in the stable region are more strin- 
gent than those obtained with the entire range of data. For 
instance, for the BaBar data @ taken as input with the 
Bern phase [l(| and the spacelike value F{— 1.60 GeV 2 ), 
we find the (r 2 ) bounds over the entire energy range to be 
(0.42 - 0.45) fm 2 and (0.42 - 0.44) fm 2 respectively. The 
weighted averages of the bounds obtained with different 
experimental data sets are very close to these values, for 
both choices of the spacelike datum, F(— 1.60 GeV 2 ) and 
F(-2.45 GeV 2 ). Thus, from Table[2]we conclude that the 
weighted averages of the results in the stable region is 
fully consistent with the bounds (|2"3"]l obtained from the 
strict intersection of the ranges in the same region, while 
the average over all the experimental points leads to the 
somewhat weaker bounds 



0.41 fni 



,.2 \ 



0.45 fni 



(25) 



The range defined by the limits given in (|23[) and the more 
conservative limits (|25[) are our parametrization-free de- 
terminations based on the recent data on the modulus. 

As we mentioned, the drastic fluctuations of the bounds 
as functions of the input energy shown in Fig. Q] are a 
signal of mutual inconsistencies in the data on modulus, 
especially in the low energy region. In order to test our 
formalism on a smoother input, we have considered the 
Gounaris-Sakurai (GS) parametrization of the CLEO data 
on the modulus j36[ , shown in Fig. [6] We generated from 
this model data at 121 discrete points, using them and the 
corresponding errors as input in our formalism for deriving 
bounds on (r 2 ). The upper and lower bounds obtained in 
this way for the two phases [l6| and [l5| are respectively 
shown in Figs. [7] and [5] as functions of the energy at which 
the input modulus was used. One may see that the bounds 
depend now smoothly on the input energy. Moreover, as 
expected, the bounds are stronger when the input is from 
the low energy region, and become weaker for the input 
from higher energies. 



0.46 



0.45 



E 0.44 



Since a strict intersection of the (r 2 ) ranges upon the 
whole energy region does not exist, we may adopt a weaker 
definition of the final allowed range, by taking the weighted N >- K 
average of the upper (lower) bounds on (r 2 ) shown in Fig. o.43 
[TJ It is reasonable to use as weights in this average the 
experimental errors e n on the measured modulus F n as 
defined in (f2"2"j). Thus we calculate for both the upper and 042 
lower bounds the weighted averages 
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Fig. 5. As in Fig. H but with the input F(— 2.45 GeV 2 ). 



For comparison, we have included the results obtained 
with the data simulated from CLEO parametrization in 
Figs. [2]- 03 which show the bounds obtained with the in- 
put from the stable region. We also present in Tables [T] 
and [2] the intersections of the individual ranges and the 
weighted averages, respectively, calculated as in the case 
of the real data. From the left part of Table [1] one may see 
that the full intersection of the ranges over the whole set of 
energies is practically empty also in this case. This shows 
that the input data on the modulus, although smooth, are 
not entirely consistent with the analyticity requirements, 
if the phase in the elastic region and the spacclike val- 
ues are imposed as in our work (recall that these pieces 
of input are taken from independent sources, Roy equa- 
tions or experiment, which provide their values with great 
precision) . 




ob , i , i , i ,_ 

0.3 0.4 0.5 0.6 

t 1 ' 2 [GeV] 



0.8 0.9 



Fig . 6. Gounaris-Sakurai parametrization of the CLEO data 
[3a| on the modulus. 



The values obtained with the parametrization of the 
modulus measured by CLEO shown in Tables [1] and [2] are 
consistent with those obtained with the data on modulus 
from the recent experiments. This exercise shows that our 



B. Ananthanarayan et al.: Shape parameters of the pion form factor 



0.45 
0.44 




Bern phase [16] + F(-1 .60 GeV 2 ) [4,5] 
Bern phase [16] + F(-2.45 GeV 2 ) [4,5] 


I 1 1 
^^^^^^^ 








0.43 








0.42 








0.41 








0.4 






^^^^^^ 


0.39 




I , I , I 


\ 

+ 
\ 

1,1, 



I I I I I I I I I I I I I 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 



t 1 ' 2 [GeV] 

Fig. 7. Upper and lower bounds on (r^) in terms of the energy 
of the input, using data simulated from the parametrization of 
CLEO data on the modulus [3f|. 
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Fig. 8. As in Fig. with the phase from [1^ . 

parametrization free method works indeed efficiently. It 
does not require an a prior parametrization of the data on 
the modulus, thereby avoiding any bias and allowing one 
to detect in a sensitive way the smallest inconsistencies in 
the input data. 



4.2 c-d analysis 

We now turn to the implications of the new data on the 
higher shape parameters c and d of the pion form factor 
at t = 0. In our previous work [2(| we have obtained an 
allowed domain in the c-d plane using as input the radius 
(r 2 ) taken from the literature, along with the phase in 
the elastic region from Fermi- Watson theorem and Roy 
equations for irir scattering and the experimental space- 
like data given in (j2"Tj) . The measurements of the modulus 
in the timelike region below the inelastic threshold were 
not used in that study. In the present work, we update 
the analysis by including these data as well, in order to 
see their impact on the determination of the higher order 
shape parameters at t = 0. 



In [20l ] we have taken as input a rather narrow range 
(r^) = (0.435 ± 0.005) fm 2 . Now we use the analysis of 
the previous section, which favours the admissible range 
(0.42 — 0.44) fm 2 , yielded by the intersection of the indi- 
vidual ranges from the stable region and by the weighted 
average of the bounds from the same region. We consider 
also the more conservative range (0.41—0.45) fm 2 obtained 
from the weighted average of the bounds from all energies 
in the elastic region. 

The allowed domain in the c-d plane is obtained from 
the inequality (|13|) written in the particular case K — 3 
and N = 2, using as input the normalization (|20[) and the 
first derivative (r 2 ) at t = 0, one spacelike value (|2ip and 
the modulus |-F(i n )| at a fixed energy t n . For a fixed input 
the allowed domain is the interior of an ellipse in the c-d 
plane. In practice, the boundary of the ellipse is given by 
the upper and lower bounds on d at each fixed c, found 
by solving a quadratic equation. The input quantities (the 
phase, the radius, the spacelike datum and the modulus 
|.F(t n )| at a fixed energy) were then varied within their 
error bars and the weakest bounds on d at each fixed c 
were taken, as explained at the end of Section [31 Finally, 
the input energy t n was varied and the intersection of the 
particular allowed domains obtained at various energies 
was obtained, by taking the smallest upper bound and 
the largest lower bound on d at each admissible c. The 
boundary of the resulting allowed domain is no longer an 
ellipse in the c-d plane, but has a more complicated shape 
found numerically. 

As a first step, we repeat the analysis done in (20j 
without input from the timelike modulus data below t- m , 
using now a more conservative range for the charge radius 
(r 2 ). We use also a different weight function p(t) in the 
condition ©: while in [2(| we adopted a special weight 
relevant for the calculation of the muon's (g — 2), we now 
make calculations with p = 1/t and p = l/y/t. The results 
are presented in Figs. [H]-[T21 separately for the narrow and 
conservative range of (r 2 ), the two input phases, Bern [l6| 
and Madrid , the two spacelike values (f2~Tj) and the two 
choices of the weight. A detailed examination shows that 
the domains obtained with various inputs are consistent 
between them. The results are also consistent with the 
ranges of c and d quoted in [20], being slightly weaker 
due to the more conservative input for the charge radius 
adopted now. 

Next, we impose also the timclikc modulus data be- 
low the inelastic threshold as an additional constraint. We 
first illustrate the impact of these data on the c-d domain 
by presenting results obtained with a specific input, then 
make a systematic investigation by exploiting the full set 
of data at all energies. 

The four panels of Fig. ll3l show the allowed domains in 
the c-d plane obtained using as input the modulus \F(t)\ 
measured by BaBar experiment [6| at four particular en- 
ergies. We worked with the weight p(t) = 1/t in ([3]). The 
Madrid phase [l5|, the first spacelike point from (f2"Tjl and 
the value (r 2 ) = 0.43 fm were used as input in the cal- 
culations. The three ellipses in each panel are obtained 
with the central value and the extreme values of the in- 
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Fig. 9. Allowed domain in the c-d plane with no input on the 
modulus below t- ln , for (t- 2 ) in the range (0.41 - 0.45) fm 2 and 
one spacelike datum from (|2ip . The weight function in @ is 
p(t) = l/t. 
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Fig. 10. As in Fig. [9] with (r 2 ) in the range (0.42 - 0.44) fm 2 



ut modulus at the ends of the error intervals quoted in 
The comparison with Figs. l9l-[T2l shows the reduction 
of the allowed domain by the additional constraint on the 
modulus. The figure illustrates also the effect of the exper- 
imental errors on the shape parameters: the three ellipses 
are quite distinct for the input at 0.701 GeV, but practi- 
cally coincide for the input at 0.901 GeV, where the errors 
are smaller. We obtained similar ellipses also for the other 
data sets, KLOE, CMD-2 and Belle. Of course, the allowed 
domain for the input modulus from a fixed energy will be 
obtained by the union of the individual ellipses obtained 
by varying the input quantities inside the error intervals, 
as explained above. 

In order to include the constraint from the modulus 
measurements at all energies, we finally have to take the 
intersection of the individual allowed domains obtained 
using as input the modulus measured at a fixed energy. 
If the measurements at different energies are consistent 
with each other, the intersection will be non empty and 
not very different from the individual domains. This is not 
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Fig. 11. As in Fig. [9] with the weight function pit) = l/\fi 
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Fig. 12. As in Fig.[TTl with (r 2 ) in the range (0.42-0.44) fm 2 



the case if we consider all the points measured by BaBar 
below the inelastic threshold. Therefore, as in the previous 
section we restrict the energy interval to the stable region 
(0.65, 0.70) GeV, for which a nonzero intersection exists. 
In Figs. [T4l and [T5l we present the analysis for two choices 
of the allowed interval for (r 2 ), the conservative range 
(0.41 - 0.45) fm 2 and the narrow range (0.42 - 0.44) fm 2 , 
respectively. The comparison with the domains shown in 
Figs. l9l- H2l shows the improvement brought by the input 
modulus below the inelastic threshold. From Figs. [HI and 
1151 we extract the conservative ranges for the higher shape 
parameters: 



(26) 



c 6 (3.86, 4.05) GeV" 4 , 
tie (10.32, 10.58) GeV" 6 , 

with a strong correlation between them. 

5 Discussion and conclusions 



We will now summarize the results of our analysis and 
present some conclusion and discussions. The phenomeno- 
logical knowledge of the pion electromagnetic form factor 
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Fig. 13. Ellipses in the c-d plane obtained with input on the 
modulus measured at several energies below t- m by BaBar 

experiment. The remaining input is described in the text. 
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Fig. 14. Allowed domain in the c-d plane using the modu- 
lus data from [(| in the range (0.65 - 0.70) GeV and (r 2 ) in 
the range (0.41 — 0.45) fm 2 . The input phase and the spacelike 
datum are specified in each panel. 



improved considerably in recent years. The information 
comes from different sources: the phase below 0.917 GeV is 
known via Fermi- Watson theorem from the P-wave phase 
shift of 7T7T scattering calculated with precision from the 
Roy equations, the modulus has been measured recently 
by high statistics experiments on e + e~ annihilation or 
T-decays, and measurements on the spacelike axis from 
electroproduction data with improved accuracy were also 
reported. 

If the modulus and the phase are both known along 
an energy interval the function can be reconstructed in 
principle everywhere, by the uniqueness of analytic con- 
tinuation. However, since the input is known with some 
uncertainty, the reconstruction is practically impossible. 
Indeed, analytic continuation is unique, but it is also an 
unstable (ill-posed) problem in the Hadamard sense fl7j ]. 
which means that the small uncertainties along the origi- 
nal range may be amplified arbitrarily much in the com- 
plex plane outside that range. Therefore, suitable tools 
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Fig. 15. As in Fig-liH with (r 2 ) in the range (0.42-0.44) fm 2 



based on analyticity must be devised for exploiting the 
available knowledge in order to make model-independent 
predictions in regions not directly accessible to experi- 
ment. 

In our analysis we applied a formalism that exploits in 
an optimal way the information on the phase and modulus 
on the unitarity cut, leading to results independent on the 
unknown phase above the inelastic threshold. Combined 
with the well-known analytic interpolation theory, it al- 
lows us to include also information or derive predictions 
on the values at points inside the analyticity domain. The 
technique has been already applied to the pion form factor 
in Refs. 0, M, H3 • 

In the present work we studied the impact of the time- 
like modulus data below the inelastic threshold on the 
determination of the charge radius (rj) and the higher 
shape parameters c and d in the Taylor expansion ([1]). 
The formalism also acts as a sensitive devise for testing 
the consistency of the various experimental data sets with 
analyticity and among themselves. To our knowledge, this 
is the first application of the high- statistics modulus mea- 
surements to the extraction of the pion charge radius and 
the higher shape parameters. 

Using as input the conditions ©, ©, ([201), (ED) and 
([%!?]) ■ we calculated from the inequality (|13p upper and 
lower bounds on the first derivative appearing in the ex- 
pansion ([T]). For a given input, the bounds are optimal 
and imply no loss of information, being obtained by solv- 
ing exactly a functional interpolation problem. Moreover, 
as already emphasized, the knowledge of the phase of the 
form factor in the inelastic region is not required. The un- 
certainties of the phase in the elastic region, the spacelike 
value and the modulus were taken into account by varying 
the input inside the error intervals and taking the largest 
allowed range. The final prediction for (rj) is given by 
the intersection of all the ranges obtained with an input 
modulus at fixed energy below the inelastic threshold. We 
have performed this analysis separately for all data sets of 
BaBar, KLOE, CMD-2 and Belle experiments. 

The analysis presented in Sec. 0] shows that, for all 
the experimental data sets, this intersection turns to be 
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empty. This shows that the data on the modulus are not 
all consistent with each other. Especially the low energy 
data lead to results that strongly fluctuate. Moreover, for 
some of the data real bounds on (r 2 ) cannot be obtained 
from (|13[) , which means that these data are in conflict with 
analyticity and the input phase. However, if we restrict the 
input to measurements of the modulus in the energy range 
(0.65 - 0.70) GeV, the results from all the experiments, 
BaBar, KLOE, CMD-2 and Belle, are consistent and lead 
to the prediction (|23|) . which we write as 

(r 2 ) g (0.42, 0.44) fm 2 . (27) 

For the higher shape parameters c and d, the ranges ob- 
tained with the BaBar data are given in ((261) . 

We recall that in [22j the mathematical formalism was 
applied in a somewhat opposite way, i.e. we adopted as 
input the value of the radius and derived bounds on the 
modulus below the cott threshold. In fact, we have used 
as input precisely the range (|27[). which we now obtained 
from the information on the modulus, and found that the 
derived bounds are consistent with the experimental data 
at higher energies, but in some disagreement at low en- 
ergies. Thus, the two analyses are perfectly consistent. In 
particular, as already mentioned, the bounds on the mod- 
ulus at low energies, derived with the input range (|27p 
for the radius, are more precise that the data. Therefore, 
they can be used for improving the evaluation of the two- 
pion contribution to the muon's magnetic anomaly. This 
problem will be studied in a future work. 

Note that we avoid the presentation of our results 
in terms of a central value and an error, for instance 
(rl) = (0.43 ± 0.01) fm 2 instead of fl23- The reason is 
that, although the probability of the parameters to be in 
the predicted allowed ranges follows from the probabili- 
ties of the error intervals inside which we varied the input 
quantities, we cannot actually attach precise probabilities 
to the specific values inside the allowed domains. Indeed, 
even for the central values of the input we obtain a range 
for the parameters of interest, not a definite prediction. 
As a statistical interpretation is difficult to give, we may 
consider that all the values in the predicted ranges are 
equally probable. 

The range (f2"T)l is consistent with most of the results 
based on specific parametrizations reported in the litera- 
ture, while for the higher shape parameters our results 
(|2"6"|) are more precise than the previous ones. We em- 
phasize again that our predictions do not rely on specific 
parametrizations of the modulus or the form factor itself. 
Extending our previous studies H?J HH, [22[, the present 
analysis illustrates the usefulness of appropriate analytic- 
ity tools in conjunction with high accuracy data for im- 
proving the description of the pion electromagnetic form 
factor. 

Acknowledgement: IC acknowledges support from the 
Program Nucleu under Contract PN 09370102/2009, and 
from Program Idci-PCE, Contract No 121/2011. 



References 

1. S.R. Amendolia et al. [NA7 Collaboration], Nucl. 
Phys. B 277 (1986) 168. 

2. C.N. Brown, C.R. Canizares, W.E. Cooper, A.M. 
Eisner, G.J. Feldmann, C.A. Lichtenstein, L. Litt, 
W. Loccrctz, V.B. Montana and F.M. Pipkin Phys. 
Rev. D 8, 92 (1973). 

3. C.J. Bebck et al.,' Phys. Rev. D17, 1693 (1978). 

4. T. Horn et al. [Jefferson Lab F(pi)-2 Collaboration], 
Phys. Rev. Lett. 97 (2006) 192001 [nucl-ex/0607005 . 

5. G. M. Hubcr et al. [Jefferson Lab Collaboration], 
Phys. Rev. C 78 (2008) 045203 |arXiv:0809.3052l 
[nucl-ex]]. 

6. B. Aubcrt et al. [BABAR Collaboration], P hys. Rev. 
Lett. 103 (2009) 231801 |arXiv:0908.3589l [hep-ex]]. 

7. F. Ambrosino et al. [KLOE Collaboration], Phys. 
Lett. B 670 (2009) 285 |arXiv:0809.3950l [hep-ex]]. 

8. F. Ambrosino et al. [KLOE Collaboration], Phys. 
Lett. B 700 (2011) 102 |arXiv:1006.53l3l [hep-ex]]. 

9. R.R. Akhmetshin et al., [CMD-2 Collaboration], 
JETP Lett. 84 (2006) 4 13 [Pisma Z h. Eksp. Teor. 
Fiz. 84 (2006) 491] [hep- ex/ 06 10016] . 

10. R. R. Akhmetshin et al. | CMD-2 Collabora tion] , 
Phys. Lett. B 648 (2007) 28 [hep-ex/0610021| . 

11. M. Fujikawa et al. [Belle Collaboration] , Phys. Rev. 
D 78 (2008) 072006 |arXiv:0805.3773l [hen-ex]]. 

12. B. Ananthanarayan, G. Colangclo, J. Gasser 
and H. Leutwyl er, Phys. Rept. 353 (2001) 207 
[hep-ph/0005297] . 

13. G. Colangelo, J. Gasser and H. Leutwy ler, Nucl. 
Phys. B 603 (2001) 125 |hep-ph/0103088] . 

14. R. Kaminski, J. R. Pelaez and F. J. Yndurain, Phys. 
Rev. D 77 (2008) 054015 |arXiv:0710. 11501 [hep-ph]]. 

15. R. Garcia-Martin, R. Kaminski, J. R. Pelaez, J. Ruiz 
de Elvir a and F. J. Yndur ain, Phys. Rev. D 83 (2011) 
074004 |arXiv:1102.2l83l [hep-ph]]. 

16. I. Caprini, G. Colangelo and H. Leutwyler, Eur. 
Phys. J. C 72 (2012) 1860 larXiv:1111.7160l [hep- 
ph]]. 

17. S. Chilli, C. Pomponiu and I. Sabba-Stefancscu, 
Phys. Rept. 17 (1975) 133. 

18. I. Caprini, Eur . Phys. J. C 13 (2000) 471 
[hep-ph/9907227] . 

19. G. Abbas, B. Ananthanarayan, I. Caprini, I. Sen- 
titemsu Imsong and S. Ramanan, Eur. Phys. J. A 
45 (2010) 389 [arXiv: 1004.42571 [hep-ph]]. 

20. B. Ananthanarayan, I. Caprini and I. S. Imsong, 
Phys. Rev. D 83 (2011) 096002 |arXiv:1102.3299l 
[hep-ph]]. 

21. B. Ananthanarayan, I. Caprini and I.S. Imsong, 
Phys. Rev. D 85 (2012) 096006 |arXiv:1203. 53981 
[hep-ph]]. 

22. B. Ananthanarayan, I. Caprini, D. Das and 
I. S. Imsong, E ur. Phys. J. C 72 (2012) 2192 
|arXiv:1209.0379l [hep-ph]]. 

23. G. Colangelo, M. Finkcmeicr and R. Urcch, Phys. 
Rev. D 54, 4403 (1996) [arXiv:hep-ph/9604279| . 

24. J. Bijnens, G. C olangelo and P. Talave ra, JHEP 
9805, 014 (1998) [arXiv:hep-ph/9805389] . 



12 



B. Ananthanarayan et al.: Shape parameters of the pion form factor 



25. J.Bijncns and P. Talave ra, JHEP 0203, 046 (2002) 
|arXiv:hep ph/020304^ . 

26. S. Aoki et al. [JLQCD Collaboration and TWQCD 
Collaboration], Phys. Rev. D 80, 034508 (2009) 
|arXiv:0905.2465l [hep-lat]] . 

27. J. F. de Troconiz a nd F. J. Yndurain, Phy s. Rev. D 
71, 073008 (2005) |arXiv:nep-ph/0402285] . 

28. P. Masjuan, S. Peris and J. J. Sanz-Cillero, Phys. 
Rev. D 78, 074028 (2008) ar Xiv:0807.4893l [hep-ph]] . 

29. T. N. Truong, [arXiv:hep -ph/9809476l 

30. R. Baratc et al. [ALEPH Collaboration], Z. Phys. C 
76, 15 (1997). 

31. G. J. Gounaris and J. J. Sakurai, Phys. Rev. Lett. 
21, 244 (1968). 

32. N.N. Meiman, Sov. Phys. JETP. 17 (1963) 830. 

33. P. Duren, Theory of H p Spaces, Academic Press, 
New York, 1970 

34. H. Leutwyler, |hep-ph/0212324 

35. C. Hanhart, P hys. Lett. B715 (2012) 170, 
larXiv:1203.6839l [hep-ph] . 

36. S. Anderson et al. [CLEO Collaboration], Phys. Rev. 
D 61 (2000) 112002 [hep-ex/9910046]. 



